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Reverberation 

Reverberation is a property of concert halls that greatly adds to the 
enjoyment of a musical performance. Sound waves propagate directly from 
the stage to the listener, and also reflect from the floor, walls, ceiling, and 
back wall of the stage to create myriad copies of the direct sound that are 
time-delayed and reduced in intensity. In this module, learn about the 
concept of reverberation in more detail and ways to emulate reverberation 
using a digital filter structure known as a comb filter. 


This module refers to LabVIEW, a software development 
environment that features a graphical programming language. 
Please see the LabVIEW QuickStart Guide module for tutorials 
and documentation that will help you: 


¢ Apply LabVIEW to Audio Signal Processing 
* Get started with LabVIEW 


* Obtain a fully-functional evaluation edition of LabVIEW 


Introduction 


Reverberation is a property of concert halls that greatly adds to the 
enjoyment of a musical performance. The on-stage performer generates 
sound waves that propagate directly to the listener's ear. However, sound 
waves also bounce off the floor, walls, ceiling, and back wall of the stage, 
creating myriad copies of the direct sound that are time-delayed and 
reduced in intensity. 


In this module, learn about the concept of reverberation in more detail and 
ways to emulate reverberation using a digital filter structure known as a 
comb filter. 


The [link] screencast video continues the discussion by visualizing the 
sound paths in a reverberant environment. The impulse response of the 
reverberant environment is also introduced. Understanding the desired 
impulse response is the first step toward emulating reverberation with a 
digital filter. 

https://youtu.be/agU2B5XgCOM (0:48) 

https://www. youtube.com/embed/agU2B5X gCOM?Prel=0 


[video] Sound paths in a reverberant environment 
and impulse response model 


Comb Filter Theory 


The comb filter is a relatively simple digital filter structure based on a 
delay line and feedback. Watch the [link] screencast video to learn how 
you can develop the comb filter structure by considering an idealized 
version of the impulse response of a reverberant environment. 
https://youtu.be/K19eykmRaho (1:36) 
https://www.youtube.com/embed/K19eykmRaho?rel=0 


[video] Development of the comb filter structure 
as a way to implement an idealized reverberant 
impulse response 


The difference equation of the comb filter is required in order to 
implement the filter in LabVIEW. In general, a difference equation states 
the filter output y(n) as a function of the past and present input values as 
well as past output values. 


Take some time now to derive the comb filter difference equation as 
requested by the following exercise. 
Exercise: 


Problem: 


Derive the difference equation of the comb filter structure shown at the 
end of the [link] video. 


Solution: 


Comb Filter Implementation 


Once the difference equation is available, you can apply the coefficients of 
the equation to the LabVIEW "IIR" (infinite impulse response) digital filter. 
The [link] screencast video walks through the complete process you need to 
convert the comb filter difference equation into a LabVIEW VI. The 
LabVIEW MathScript node creates the coefficient vectors. Once the comb 
filter is complete, its impulse response is explored for different values of 
delay line length and feedback gain. 


Please follow along with the video and create your own version of the comb 
filter in LabVIEW. Refer to TripleDisplay to install the front-panel indicator 
used to view the signal spectrum. 

https://youtu.be/LCZ1kaP_h6s (11:16) 
https://www.youtube.com/embed/LCZ1kaP_h6s?rel=0 


[video] Implementing the comb filter in 
LabVIEW; exploration of the impulse response 
as a function of delay line length and feedback 

gain 


Loop Time and Reverb Time 


As you have learned in previous sections, the comb filter behavior is 
determined by the delay line length N and the feedback coefficient g. From 
a user's point of view, however, these two parameters are not very intuitive. 


Instead, the comb filter behavior is normally specified by loop time 7 and 
reverb time denoted 79. Reverb time may also be written as Rr¢go . Loop 
time indicates the amount of time necessary for a given sample to pass 
through the delay line, and is therefore the delay line length N times the 
sampling interval. The sampling interval is the reciprocal of sampling 
frequency, so the loop time may be expressed as in [link]: 

Equation: 


Reverb time indicates the amount of time required for the reverberant 
signal's intensity to drop by 60 dB (dB = decibels), effectively to silence. 
Recall from the [link] video that the comb filter's impulse response looks 
like a series of decaying impulses spaced by a delay of N samples; this 
impulse response is plotted in [link] with the independent axis measured in 
time rather than samples. 
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Comb filter impulse response 


Take a few minutes to derive an equation for the comb filter feedback gain 
"9" as a function of the loop time and the reverb time. The following pair of 
exercises guide you through the derivation. 

Exercise: 


Problem: 


Given the comb filter impulse response pictured in [link], derive an 
equation for the reverb time T7¢ in terms of the loop time 7 and the 
comb filter's feedback constant g. Hint: Recall the basic equation to 
express the ratio of two amplitudes in decibels, i.e., use the form with a 
factor of 20. 


Solution: 
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Exercise: 


Problem: 


Based on your previous derivation, develop an equation for the comb 
filter gain "g" in terms of the desired loop time and reverb time. 


Solution: 
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Exercise: 


Problem: 


To finish up, derive the equation for the comb filter delay "N" in terms 
of the desired loop time. 


Solution: 


N=T{¢4 


Now, return to your own comb filter VI and modify the front-panel controls 
and LabVIEW MathScript node to use loop time and reverb time as the 
primary user inputs. Experiment with your modified system to ensure that 
the spacing between impulses does indeed match the specified loop time, 
and that the impulse decay rate makes sense for the specified reverb time. 


Comb Filter Implementation for Audio Signals 


In this section, learn how to build a comb filter in LabVIEW that can 
process an audio signal, specifically, an impulse source. Follow along with 
the [link] screencast video to create your own VI. The video includes an 
audio demonstration of the finished result. As a bonus, the video also 
explains where the "comb filter" gets its name. 


[video] Building a LabVIEW VI of a comb filter 
that can process an audio signal 


Next, learn how you can replace the impulse source with an audio .wav file. 
Speech makes a good test signal, and the [link] screencast video shows how 
to modify your VI to use a .wav audio file as the signal source. The speech 
clip used as an example in the video is available here: speech.wav (audio 
courtesy of the Open Speech Repository, 
www.voiptroubleshooter.com/open_ speech; the sentences are two of the 
many phonetically balanced Harvard Sentences, an important standard for 
the speech processing community). 


[video] Modifying the LabVIEW VI of a comb 
filter to process a .wav audio signal 
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Schroeder Reverberator 

The Schroeder reverberator uses parallel comb filters followed by cascaded 
all-pass filters to produce an impulse response that closely resembles a 
physical reverberant environment. Learn how to implement the Schroeder 
reverberator block diagram as a digital filter in LabVIEW, and apply the 
filter to an audio .wav file. 


This module refers to LabVIEW, a software development 
environment that features a graphical programming language. 
Please see the LabVIEW QuickStart Guide module for tutorials 
and documentation that will help you: 


* Apply LabVIEW to Audio Signal Processing 
* Get started with LabVIEW 


* Obtain a fully-functional evaluation edition of LabVIEW 


Introduction 


Reverberation is a property of concert halls that greatly adds to the 
enjoyment of a musical performance. The on-stage performer generates 
sound waves that propagate directly to the listener's ear. However, sound 
waves also bounce off the floor, walls, ceiling, and back wall of the stage, 
creating myriad copies of the direct sound that are time-delayed and 
reduced in intensity. 


In the prerequisite module Reverberation, you learned how the comb filter 
structure can efficiently create replicas of a direct-path signal that are time 
delayed and reduced in intensity. However, the comb filter produces 
replicas that are time delayed by exactly the same amount, leading to the 
sensation of a pitched tone superimposed on the signal. Refer back to 


Reverberation to hear an audio demonstration of this effect. Put another 
way, the impulse response of the comb filter contains impulses with 
identical spacing in time, which is not realistic. 


The Schroeder reverberator (see "References" section) uses a 
combination of comb filters and all-pass filters to produce an impulse 
response that more nearly resembles the random nature of a physical 
reverberant environment. 


This module introduces you to the Schroeder reverberator and guides you 
through the implementation process in LabVIEW. As a preview of what can 
be achieved, watch the [link] screencast video to see and hear a short 
demonstration of a LabVIEW VI that implements the Schroeder 
reverberator. The speech clip used in the video is available here: 

speech. wav (audio courtesy of the Open Speech Repository, 
www.voiptroubleshooter.com/open_ speech; the sentences are two of the 
many phonetically balanced Harvard Sentences, an important standard for 
the speech processing community). 

https://youtu.be/ZgxjjBsB7_M (0:00) 
https://www.youtube.com/embed/Z¢gxjjBsB7_ MPrel=0 


[video] Demonstration of the Schroeder 
reverberator as implemented in LabVIEW 


Structure of the Schroeder Reverberator 


The [link] screencast video presents the structure of the Schroeder 
reverberator and describes the rationale for its design. 
https://youtu.be/k1SorBkeglo (0:00) 
https://www.youtube.com/embed/k1SorBkeglo?rel=0 


[video] Structure of the Schroeder 
reverberator and rationale for its design 


All-Pass Filter 


The Schroeder reverberator uses all-pass filters to increase the pulse 
density produced by the parallel comb filters. You perhaps are familiar with 
the frequency response of an all-pass filter: its magnitude response is unity 
(flat) for all frequencies, and its phase response varies with frequency. For 
example, the all-pass filter is used to create a variable fractional delay as 
described in Karplus-Strong Plucked String Algorithm with Improved Pitch 
Accuracy. 


From an intuitive standpoint it would seem that a flat magnitude response in 
the frequency domain should correspond to an impulse response (time- 
domain) containing only a single delta function (recall that the impulse 
function and a constant-valued function constitute a Fourier transform pair). 
However, the all-pass filter's impulse response is actually a large negative 
impulse followed by a series of positive decaying impulses. 


In this section, derive the transfer function and difference equation of the 
all-pass filter structure shown in [link]. Note that the structure is 
approximately of the same level of complexity as the comb filter: it 
contains a delay line of N samples and an extra gain element and summing 
element. 


All-pass filter structure 


Exercise: 


Problem: 


Determine the transfer function H(z) for the all-pass filter structure of 
[ink]. 


Solution: 
<= + —N 
H(2) = 7H 
Exercise: 
Problem: 


Based on your result for the previous exercise, write the difference 
equation for the all-pass filter. 


Solution: 


y(n) = —gz(n) + a(n — N) + gy(n — N) 


All-Pass Filter Impulse Response 


As described in the video of [link], two all-pass filters are placed in cascade 
(series) with the summed output of the parallel comb filters. An 
understanding of the all-pass filter impulse response reveals why a cascade 
connection increases the pulse density of the comb filter in such a way as to 
emulate the effect of natural reverberation. 


The [link] screencast video derives the impulse response of the all-pass 
filter; the loop time and reverb time of the all-pass filter are also presented. 
https://youtu.be/Utvvsr_dQ02c (0:00) 

https://www. youtube.com/embed/Utvvsr_d02c?rel=0 


[video] Derivation of the all-pass filter impulse 


response 


The [link] screencast video demonstrates the sound of the all-pass filter 
impulse response compared to that of the comb filter. Moreover, the audible 
effect of increasing the comb filter pulse density with an all-pass filter is 
also demonstrated in the video. 


Download the LabVIEW VI presented in the video: apfdemo.zip Refer to 
TripleDisplay to install the front-panel indicator required by the VI. 
https://youtu.be/h_bDexyaizyY (0:00) 
https://www.youtube.com/embed/h_ bDexya1zY ?rel=0 


[video] Audio demonstration of the all-pass filter 
impulse response combined with comb filter 


Project: Implement the Schroeder Reverberator in LabVIEW 


As described in the [link] screencast video, two all-pass filters are placed in 
cascade (series) with the summed output of four parallel comb filters. The 
video of [link] explains how the all-pass filter "fattens up" each comb filter 
output impulse with high density pulses that rapidly decay to zero. 
Selecting mutually-prime numbers for the loop times ensures that the comb 
filter impulses do not overlap too soon, which further increases the effect of 
randomly-spaced impulses. 


The table in [link] lists the required reverb times (T60) and loop times (tau) 
of the Schroeder reverberator. Note that the comb filters all use the same 
value (the desired overall reverb time). 


Reverb time and loop time values for the comb filters 
and all-pass filters of the Schroeder reverberator 


The [link] screencast video provides everything you need to know to build 
your own LabVIEW VI for the Schroeder reverberator. 


Download the .wav reader subVI mentioned in the video : WavRead.vi. 
https://youtu.be/nW7haGUkm8A (0:00) 
https://www. youtube.com/embed/nW7haGUkm8A ?rel=0 


[video] Suggested LabVIEW techniques to build 
your own Schroeder reverberator 
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Localization Cues 

Learn about two localization cues called interaural intensity difference (IID) 
and interaural timing difference (ITD), and learn how to create a LabVIEW 
implementation that places a virtual sound source in a stereo sound field. 


This module refers to LabVIEW, a software development 
environment that features a graphical programming language. 
Please see the LabVIEW QuickStart Guide module for tutorials 
and documentation that will help you: 


¢ Apply LabVIEW to Audio Signal Processing 
* Get started with LabVIEW 


* Obtain a fully-functional evaluation edition of LabVIEW 


Introduction 


Our enjoyment of synthesized or recorded sound is greatly enhanced when 
we perceive the actual locations of the various musicians on stage. In a high 
quality stereo recording of a jazz trio, for example, you can tell that the 
drummer is perhaps located slightly to the left of center, the saxophonist is 
on stage right and the bass player is on stage left. If you have ever flipped 
on the "mono" (monaural) switch on your stereo amplifier, then you know 
that the resulting sound is comparatively unpleasant, especially when 
wearing headphones. 


We live in a three-dimensional soundfield, and our ears continually 
experience slightly different versions of any given sound. These variations 
allow us to place (or localize) the sound source, even when we cannot see 
it. 


In this module, learn about two localization cues called interaural 
intensity difference and interaural timing difference, abbreviated IID 
and ITD, respectively. Each cue relies on presenting a slightly different 
version of a sound to each ear. 


Interaural Intensity Difference (IID) 


The interaural intensity difference (IID) localization cue relies on the fact 
that an off-center source produces a higher intensity sound in one ear 
compared to the other. This technique is also called intensity panning, and 
is most effective when the listener is wearing headphones. When a multi- 
speaker arrangement is used in a room, a given sound propagates to both 
ears. Lower frequencies have longer acoustic wavelengths and are therefore 
less directional than higher frequency sounds. Consequently, the IID effect 
works better at higher frequencies above 1400 Hz (Dodge and Jerse, 1997). 


The [link] screencast video continues the discussion by deriving the 
equations needed to implement the IID effect using a two-speaker 
arrangement. 

https://youtu.be/TOidBG7mi-w (7:48) 

https://www. youtube.com/embed/TOidBG7mi-w?rel=0 


[video] Development of the equations needed to 
implement the interaural intensity difference 
(IID) localization cue 


Interaural Timing Difference (IID) 


The interaural timing difference (ITD) localization cue relies on the fact 
that sound waves from an off-center source arrives at one ear slightly after 
the other ear. We can perceive this slight difference in time down to about 
20 microseconds (Dodge and Jerse, 1997), and this difference helps us to 
place the sound source either to the left or right. The ITD cue works best in 
the lower frequency range of 270 to 500 Hz (Dodge and Jerse, 1997). 


The [link] screencast video continues the discussion by deriving the 
equation needed to implement the ITD effect using a two-speaker 
arrangement. 

https://youtu.be/tOP2UQP2U0k (2:34) 

https://www. youtube.com/embed/t0P2uQP2U0k?rel=0 


[video] Development of the equation needed to 
implement the interaural timing difference 
(ITD) localization cue 


Project: Implement the IID and ITD Localization Cues in 
LabVIEW 


You can easily experiment with both the IID and ITD localization cues in 
LabVIEW. The cues are probably easier to perceive when you choose a 
speech signal for your source. 


The ITD cue requires a delay or time shift between the two stereo channels. 
The delay line can be constructed as a digital filter with the transfer 
function shown in [link]: 

Equation: 


Stated as a difference equation, the filter is y(n) = a(n — N), which states 
that the output is the same as the input but delayed by N samples. The 
forward (or "b") coefficients are therefore zero for all but by = 1, and the 
reverse (or "a") coefficients are zero for all but ay = 1. 


The [link] screencast video provides LabVIEW coding tips to implement 
the equations and to generate a stereo audio signal. 


https://youtu.be/Ma4aNsfg15s (7:25) 
https://www. youtube.com/embed/Ma4aNsfq15s?rel=0 


[video] LabVIEW coding tips for the ITD and 
ITD localization cues 
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